Upper limits for the photoproduction cross section for the $ (1860) pentaquark state 

off the deuteron 
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We searched for the "I> (1860) pentaquark in the photoproduction process off the deuteron in 
the H~7r~ decay channel using CLAS. The invariant mass spectrum of the H~7r~ system does not 
indicate any statistically significant enhancement near the reported mass M = 1.860 GeV. The 
statistical analysis of the sideband-subtracted mass spectrum yields a 90% confidence level upper 
limit of 0.7 nb for the photoproduction cross section of $ (1860) with a consecutive decay into 
H~7r~ in the photon energy range 4.5 GeV < < 5.5 GeV. 

PACS numbers; 14.20.Jn, 14.20.Pt, 13.60.Rj, 12.39.-x 



INTRODUCTION 

Narrow bound states of four quarks and one anti-quark 
have been the focus of intense searches since the report 
by the LEPS collaboration of a positively charged baryon 
called the 8+, with S — +1 and a mass of 1.54 GeV m . 
This "exotic" combination of quantum numbers cannot 
be accommodated within the simple quark model, which 
assumes that all baryons are built out of three quarks. 
Exotic states of this kind have been predicted within 
the Chiral Soliton Model 0| as part of a spin 1/2 anti- 
decuplet of baryons. The anti-decuplet of "pentaquarks" 
contains three explicitly exotic states, whose quantum 
numbers require a minimal quark content of four quarks 
and one anti-quark. Reference describes the exper- 
imental situation for 8+ (1540) searches, and a concise 
summary of the current state of pentaquarks can also be 
found in the Particle Data Group (PDG) review |4|. 

The two other exotic states of the anti-decuplet have 
charge Q — —2 (quark content of ddssu) and Q = +1 
(quark content of uussd). Their strangeness is 5 = —2, 
but they have isospin 3/2, in contrast to normal cas- 
cade states with isospin 1/2. The Particle Data Group 
Q has assigned the name of $(1860) to the four states 
in the strangeness S — —2 sector of the anti-decuplet. 
The NA49 collaboration has reported evidence for the 
strangeness S = —2 pentaquark $ and the $° at 
a mass of 1.862 GeV [5|]. This measurement was con- 
ducted in p+p collisions at a center-of-mass energy -^s — 
17.2 GeV, and the states were reconstructed from their 
decays into the ground state cascades, $ — >■ S^tt^ 
and -> S^7r+. We also note that many experiments 
, some of which represent a much larger statistical 
sample^have not been able to confirm the NA49 obser- 
vation [5[. 

Guidance for where and how to search for cascade pen- 
taquarks is very sparse. The mass scale for the 8+ can 
be estimated to be about the mass of the nucleon (0.94 
GeV) plus the mass of the kaon (0.5 GeV). The mass 
of the cascade pentaquark contains an additional strange 



Source 


Mass (GeV) 


Chiral-Soliton Model [2, 20] 
Chiral- Soliton Model [23] 
Diquark Model [21, 24] 
Diquark-Triquark Model [22, 25] 
Experiment NA49 [5] 


2.07 [1.86] 
1.79-1.97 
1.75 
1.783 
1.862 ± 0.002 



TABLE I: Selected representative model expectations for the 
mass of the cascade pentaquark. Note that the first entry 
gives the initial prediction of the Chiral-Soliton model of 2.07 
GeV, followed in square brackets with the adjusted model to 
the experimental value of NA49. 



quark, which naively would lead to 1.89 GeV, assum- 
ing that the strange and anti-strange quarks have the 
same mass of about 0.45 GeV. The quark model predic- 
tions vary depending on the amount of mixing between 
the anti-decuplet and octet members, as well as the es- 
timated size of the color-spin hyperfine interaction be- 
tween quarks. We also note that some of the models 
[20I - I22I have used the experimental reports of either the 
e+(1540) or $(1860), or both, to set the scale, thus they 
are not entirely unbiased. Table H] summarizes represen- 
tative model predictions for the masses of the cascade 
pentaquark that range from 1.75 to 2.07 GeV. 

This experiment was mounted at Jefferson Lab to 
search for the first time for the $ exotic pentaquark 
state in real photoproduction off a neutron target with 
the subsequent decay into a final state containing three 
pious and one proton. Early experimental reports on 
the 8^ suggested that photon beams were a rich source 
of pentaquarks, and one calculation predicted that the 
production of the $ off the neutron was an order of 

\2l\ . Liu and 



magnitude larger than off the proton [26 
collaborators computed the photoproduction cross sec- 
tion a{-fn -> K+K+^ ) and a{-fp -J> K^K'^^'^), which 
correspond to similar reaction channels for the neutron 
and proton. At = 5 GeV the estimated cross sec- 
tion (7(771 — > K^K^^ ), assuming positive parity for 
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the exotic states, is between 0.4 and 1.5 nb, depending on 
the value of the gK' ne couphng [26| . There exists a large 
range of predictions for the decay widths and branching 
ratios of these exotic states [13, [H, [2^ , but the domi- 
nant decay mode is expected to lead to the ground state 
cascade — > tt^S^. In addition, the bias is that the 
states are very narrow and therefore long-lived, which 
is particularly interesting, because they are above the 
free particle decay thresholds. Therefore, our search was 
targeted to identify states with intrinsic widths that are 
smaller than the experimental resolution. 



EXPERIMENT 



The purpose of this experiment is to search specifically 
for the $ (1860) state of the spin-i anti-decuplet with 
the CLAS detector [s^ in Hall B at Jefferson Lab in a 
photoproduction experiment. The acceptance and res- 
olution of CLAS is better for charged than for neutral 
particles. The most promising topology for our experi- 
ment results from the decay sequence 



-(tt-A) 



"tt (tt p) . (1) 



The bremsstrahlung photon beam produced by a 
5.77 GeV electron beam interacts with the deuteron tar- 
get, producing a large variety of final states. The outgo- 
ing particles are detected and reconstructed in the CLAS 
detector. The energy and the interaction time of the ini- 
tial photon is determined by registering the electron in 
the Hall B photon tagging facility ^31] . 

The analysis strategy is to directly reconstruct the de- 
cay sequence ([T|) from the final state particles detected 
in CLAS. First, we identify the A(1116) using the proton 
and a tt^. Then we search for the S~(1321) by combin- 
ing the A(1116) with another negative pion. Finally, we 
analyze the invariant mass of the S~(1321)7r~ composite 
system to search for the $ pentaquark state. 

This CLAS experiment collected data during 40 calen- 
dar days at the end of 2004 and the beginning of 2005, 
amounting to approximately 25 pb~^ integrated luminos- 
ity in the tagged photon energy range 4.5 GeV < Ej < 
5.5 GeV. In order to achieve an adequate experimental 
sensitivity in a reasonable amount of time, we operated 
at an instantaneous photon flux significantly larger than 
ever used before with CLAS. The experimental data was 
carefully analyzed and cross-checked against known cross 
sections to properly take rate effects into account. 



APPARATUS 

Hall B at Jefferson Lab houses a photon-tagging sys- 
tem [sT] to conduct experiments with real photons. This 
facility allows for absolute cross section measurements 




FIG. 1: Three dimensional view of CLAS showing the 
three regions of Drift Chambers (Rl -R3), Cherenkov Coun- 
ters(CC), the Time-Of-Flight system (TOF) and the Elec- 
tromagnetic calorimeter (EC) (see text for details). On this 
picture, the photon beam travels from the upper-left corner 
to the lower-right corner. 



over a broad energy range of the incoming photons. The 
bremsstrahlung photon beam is produced by the electro- 
magnetic radiation of the primary electron beam in a thin 
(^ 5 X 10^^ r.l.) radiator. For this experiment, we used 
the tagged bremsstrahlung beam in Hall B incident on a 
40 cm long and 4 cm diameter liquid-deuterium target, 
which was located on the beam axis 50 cm upstream of 
the center of the CLAS detector. 

CLAS (see Fig. [ij is a nearly 47r detector that is well- 
suited to study reactions into final states with multiple 
charged particles. The magnetic field of CLAS [S^l is 
provided by six superconducting coils, which produce 
an approximately toroidal field in the azimuthal direc- 
tion around the beam axis. The regions between the 
cryostats are instrumented with six identical detector 
packages, also referred to as "sectors" . Each sector con- 
sists of four Start Counter (ST) paddles fs^] mainly used 
for triggering purposes, three regions of Drift Chambers 
(Rl, R2, and R3) [s^ to determine the trajectories of 
the charged particles, Cerenkov Counters (CC) [sl] for 
electron identification. Scintillator Counters (SC) [HI 
for charged particle identification based on the Time-Of- 
Flight (TOF) method, and Electromagnetic Calorimeters 
(EC) [36] used for electron identification and detection of 
neutral particles. 

The CLAS detector provides a ^ ~ 0.6% momentum 
resolution [sO] and up to 80% of 47r solid-angle cover- 
age. The efficiency for detection and reconstruction of 
charged particles in fiducial regions of CLAS is greater 
than 95%. The combined information from the track- 
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FIG. 2: [Color online] Mass spectrum of the p-K^ system. 
The yellow shaded area indicates the mass range used in this 
analysis. The red curve is the fit to a Gaussian peak with a 
polynomial background. The dashed vertical line shows the 
PDG 51] value for the mass of A(1116). 



ing in Drift Chambers and Scintillator Counters allows 
us to reliably separate protons from positive pions for 
momenta up to 3 GeV. 

ANALYSIS 
Event selection 

One of the main goals of the analysis procedure is to 
select events corresponding to the reaction 

-/d^<i>—X, (2) 

where we consider the decay sequence ([1]), and kaons in 
the final state are not required to be reconstructed. 

The A(1116) candidates are identified by considering 
every pair of positive and negative tracks with a hypothe- 
sis that these are the proton and the negative pion from a 
A(1116) decay using timing information from the scintil- 
lation counters and momentum and vertex information 
from tracking. To select the A-hyperon, an invariant 
mass cut 1.1108 GeV < Mp„- < 1.1202 GeV is applied 
as shown in Fig. [51 Because the decay products originate 
from the same point in space, a 5 cm cut is applied on 
the Distance-Of-Closest-Approach (DOCA) for the two 
tracks. We define DOCA as the length of the shortest line 
segment connecting the trajectories of these two tracks 
in the vicinity of the CLAS target. The detector reso- 
lution for the DOCA between the proton and the tt~ is 

1.5 cm. This cut reduces the contributions from pir" 
pairs that do not come from the A(1116) decay. If there 
is more than one A(1116) candidate, we choose the best 
pair based on the combined information from the match- 
ing of the invariant mass and the DOCA between the two 
tracks. 




0.5 1 1.5 2 2.5 3 3.5 4 



Ph(GeV) 

FIG. 3: [Color online] /3 versus momentum of the negative 
hadrons. The pions from decays of the A candidates are not 
included in this plot. The magenta lines indicate the vr" iden- 
tification cuts applied in this analysis. 



After selecting the best candidate pair for the A(1116), 
we proceed with combining it with the remaining nega- 
tive pions in the event, which are identified using time- 
of-flight and tracking information, the event start time 
determined from the vertex time of the already recon- 
structed A(1116), and the reference RF-time from the 
accelerator's injector. Figure [3] illustrates the negative 
pion identification used in this analysis. The main band 
corresponds to the negative pions. The magenta lines 
show the cuts applied to select the remaining 7r~'s in the 
event. For further analysis, we require that an event con- 
tains at least two more negative pions in addition to the 
TT^ in the A(1116) pair. Therefore, one can have multi- 
ple combinations of (Att" ) pairings, and we considered all 
combinations of the A and each of the remaining pions in 
the event whose track separation in space (DOCA) was 
sufhciently small. This treatment of the Att" pairs may 
lead to multiple entries in the background, but counts 
the correct pairing in the cascade peak only once. 

We also require that the time reconstructed from the 
event in CLAS matches the time of the interaction de- 
termined by the information from the photon tagger. If 
there is at least one photon detected in the tagger that 
can provide enough energy for the d{'y,'E,~'K~) K'^ K~^p 
reaction for the measured kinematics, and it is registered 
within ±3.2 ns of the interaction time from CLAS, then 
this event is kept in the data sample. Table |TT] summa- 
rizes the event selection cuts used in this analysis. 

Figure 3] shows the invariant mass spectrum of the 
Att" pairs after the cuts from Table |TT1 with the ex- 
ception of the cascade mass cut and the DOCA cut in 
the last two rows. The red dashed and blue dotted 
lines show the positions of the nominal mass for the 
S-(1321) and I]-(1385), respectively Q. We apply the 
1.3175 GeV < A/att- < 1-3265 GeV cut illustrated by 
the shaded area to select events with cascade hyperons. 
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DOCA(A) 


DOCA(A) < 5.0 cm 


Mass(A) 


1.1108 GeV < Mk < 1.1202 GeV 


DOCA(H) 


DOCA(H) < 4.5 cm 


/3 vs p 


p-dependent cut shown in Fig. [3] 


Tagger time 


-3.2 ns < - T^tx < +3.2 ns 


Missing Mass 


MA'U-^- > 2Mk + Mp 


Mass(H) 


1.3175 GeV < Ma^- < 1.3265 GeV 


DOCA(H7r) 


D0CA(H7r) < 4.5 cm 



TABLE II: Summary of the event selection. 




Ma ..(GeV) 

FIG. 4: [Color online] Invariant mass of Avr" pairs. The 
red dashed and blue dotted lines mark the positions of the 
nominal mass for =-(1321) and E-(1385), respectively [J. 
The shaded area shows the mass range used in this analysis. 



Figure [5] shows the invariant mass of Att-tt- after the 
cascade mass cut and the cut on the DOCA between the 
candidate and the negative pion described in Table llll 
There is no statistically significant structure near the re- 
ported mass of the (1860). We use the sideband sub- 
traction method to account for the background contribu- 
tion in the Att-tt- mass spectrum coming from events 
under the cascade peak in Fig. SI The Att-tt" spectra 
from the mass ranges 1.300 GeV < A/a^- < 1.310 GeV 
and 1.335 GeV < Ma^- < 1.345 GeV normahzed to the 
number of background events under the cascade peak are 
subtracted from the mass spectrum in Fig. [5] to obtain 
the sideband-subtracted 'E~tt~ mass distribution shown 
in Fig. m 

Detector simulation 

In order to relate the experimental yields to cross sec- 
tions, acceptance correction factors were calculated using 
the Monte-Carlo method. The GEANT-based detector 
simulation package incorporates the survey geometry of 
CLAS, realistic response of drift chamber and scintilla- 
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M^^-r (GeV) 

FIG. 5: [Color online] Invariant mass of Att-tt- after applying 
the cascade mass cut. The dashed line marks the position of 
the reported mass of $ (1860). 
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M,.,.(GeV) 

FIG. 6: [Color online] Number of H-tt" events per 6 MeV 
mass bin. The error bars indicate the statistical uncertainties. 
The dashed curve shows the fit to a Gaussian peak above a 
polynomial background. The center of the Gaussian in this 
plot is fixed at the center of the bin at 1.813 GeV. 



tion counters, as well as documented inefficiencies due 
to dead wires and malfunctioning photomultiplier tubes. 
Because CLAS is a complex detector covering almost a 
Att solid angle, it is virtually impossible to separate the 
efficiency calculations from the geometrical acceptance 
calculations. In this paper, the term acceptance correc- 
tion refers to a combined correction factor due to the 
geometry of the detector and the inefficiencies of the de- 
tection and reconstruction. It is defined as the ratio of 
the number of reconstructed Monte-Carlo events to the 
number of simulated events in each given kinematic bin. 

The event sample used in the acceptance calculation 
was generated using a phase-space generator with an 
event configuration in the first row in Table IIIII with- 
out any physics background. Figure [7] shows the invari- 
ant mass of the Att^ system from the data (a) and from 
the simulations (b). The simulated event sample does 
not contain any background because we are interested 
only in determining the acceptance and efficiency for the 
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FIG. 7: [Color online] Mass spectrum of Avr" from data (a) 
using restrictive cuts (see text for details), and Avr" mass 
from the GEANT-based simulations (b). 



Row 
# 


Production 
Model 


Relative 
Acceptance 


1 




1.00 


2 


^ K+K+$--ps K+K+E'Tv-ps 


+ 1.24 


3 


7d^iC+E-(2650)ps 

KfK+$--ps KjK+E-Tv-ps 


+1.47 


4 


-yd K+K+^--ps K+K+E-TT-ps 


+ 1.07 



TABLE III: Event configurations used in the Monte-Carlo 
generator for the model dependence studies. The elementary 
production is assumed to be from a neutron in a deuterium 
target and the proton, ps, is treated as a spectator in the 
reaction. 



events which contain cascade hyperons. The acceptance 
and efficiency corrections calculated for events containing 
S~7r~ are applied to the sideband-subtracted spectrum. 
Event selection criteria for the plots in Fig. [7] are stricter 
than the nominal cuts in this analysis to enhance the 
signal-to-background ratio for the ground state cascade 
peak for a visual comparison of the data with the simu- 
lations. 



Model dependence 

Because we do not know how the cross section of the 
<I> (1860) photoproduction depends on kinematics, and 
since we integrate over all of the kinematic variables, 
our estimate of the CLAS acceptance depends on the 
choice of the distribution of the events over the accessi- 
ble phase space. In order to estimate the uncertainty of 
the acceptance due to the model dependence, we stud- 
ied four event configurations; the relative acceptances 
are given in Table IIIII The events were generated us- 
ing a software package that includes essentially no spe- 
cific dynamics and mainly simulates events according to 




FIG. 8: [Color online] Acceptance of CLAS versus the invari- 
ant mass of the E~tt~ system, determined from the 'yd — >■ 
K'^E~-K~Ps phase space simulation, using a GEANT- 
based detector simulation package GSIM. The red line indi- 
cates the position of the peak seen by the NA49 collaboration. 



phase-space probabilities. The reaction in row 1 of Ta- 
ble IIIII is a four-body phase-space uniform distribution 
complemented with a spectator proton with Fermi mo- 
mentum smearing according to Ref. [37]. For the pro- 
cess in row 2, a hypothetical I]~(2650) was implemented 
in the event generator, which decays into 7^+$ (1860) 
with a total width of 35 MeV. The $ — (1860) is simu- 
lated as a particle with an infinitely narrow width at mass 
M$ = 1.862 GeV that decays only through the S^tt^ 
channel. The reactions in rows 2 and 4 are simulated 
according to two- and three-body phase space, respec- 
tively, with a Fermi-smeared spectator proton momen- 
tum spectrum according to Ref. JiTj . The process in row 
3 of Table IIIII is simulated according to two-body phase- 
space but with an additional exponential ^-dependence 
for the with a i-slope parameter b = 2.6 GeV~^. 
The flat phase space provides an estimate of the accep- 
tance for s-channel processes, while including a steeper 
t-dependence allows us to consider the processes going 
through the t-channel exchanges. For the processes in 
rows 2 to 4, the acceptance is estimated at a fixed mass 
M=;-^- — 1.862 GeV, since the process explicitly in- 
cludes the $^^(1860) decay. Using this table, we cal- 
culate the root-mean-square (RMS) of the differences in 
the acceptance at AI^-^- = 1.862 GeV, and we assign a 
relative uncertainty of ^21% due to model depen- 
dence. This is the largest source of uncertainty in the 
determination of the cross sections and its upper limits. 

The acceptance of CLAS versus the invariant mass of 
the S^TT^ system calculated using the process in row 1 of 
Table lllll is shown in Fig. |8l At the expected $ mass 
M — 1.862 GeV, the acceptance does not exhibit any spe- 
cial features and is ^ 0.3%. The artificial enhancement 
near mass M = 1.53 GeV is due to pion combinatorics 
and the mass cut to select S~(1321) events. 
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Trigger conditions and normalization 

The process of interest for this experiment is — > 
p7r~7r~7r~X, where three charged pions and a proton are 
detected in CLAS. To be able to run at higher photon- 
deuteron luminosities, we used a highly selective trigger, 
so that the data acquisition system could cope with the 
data rate. After extensive studies, we decided to require 
at least three charged particles in three different sectors 
of CLAS to be detected in the Start Counter and TOF 
system in the main trigger. But we also took data with a 
two-sector trigger to quantify trigger inefficiencies, which 
was prescaled by a factor of 5 to 20, depending on the 
running conditions. The trigger required a time coinci- 
dence of tracks in CLAS with a signal from the photon 
tagger, signaling the production of a photon in the energy 
range of 4.5 to 5.5 GeV. 

The primary trigger condition suffered from inefficien- 
cies at high luminosity. These inefficiencies were deter- 
mined empirically by studying the luminosity dependence 
of A production and also by comparing these yields with 
those measured with the two-sector trigger, corrected for 
the prescale factor. The study used events containing an 
identified A-hyperon, decaying to a proton and a pion, 
and two additional reconstructed tracks, similar to our 
sample of signal events. The average inefficiency of the 
three-track trigger varied linearly with the electron beam 
current up to 35%, with an average of 25% at the nominal 
current of 30 nA. The yields normalized by this factor 
correspond to those obtained from the two-track trigger, 
which did not exhibit any dependence on luminosity. 

The overall normalization of the experiment was 
checked using the c?(7, 7r~ A++)n reaction with the de- 
tection of the 7r+7r~p final state. The cross section of 
this process is expected to be mostly dominated by pho- 
toproduction off a quasi-free proton 7p — tt^A"*""^. Al- 
though final state interactions (FSI) contribute to this 
process, we do not expect their impact to be signifi- 
cant within the precision required for this purpose. Us- 
ing the same analysis and assumptions, the measured 
cross sections for d(7, tt" A++)n obtained from the cur- 
rent data was compared to the p(7, tt" A"*""*") cross section 
from a different CLAS run period where these trigger in- 
efficiencies were not present. The cross sections from 
the two CLAS data sets for the photon energy range of 
4.5 GeV < < 5.0 GeV differed by 9 - 14%. The 
agreement of our measurements with the published data 
on 7p ^ 7r-A++ from SAPHIR at E^=2.b GeV (13 is 
better than 10%, which indeed indicates that FSI con- 
tributions are negligible. A more detailed analysis of 
the (j(7, ■K~^TT~p)n reaction and its cross section is under 
way [39|. Based on these comparisons, we have assigned 
the normalization uncertainty of ±15% for the presented 
data. 



Source of Uncertainty Uncertainty (a) 

Model dependence of acceptance 21% 
Flux and trigger efficiency 15% 
Total in quadrature 26% 



TABLE IV: Sources and the values of the relative systematic 
uncertainties. The overall systematic uncertainty calculated 
as the square root of the quadrature sum is ±26%. 

Systematic uncertainties 

The final invariant mass spectrum of S^tt^ was stud- 
ied for various values of the selection parameters for the 
A, S, and 7r~. These variations did not result in any 
qualitative change of the invariant mass distributions of 
Att^tt^ or S~7r~. We checked the sensitivity of the re- 
sults with respect to the following parameters: selection 
parameters for the A-candidates, cuts on time matching 
between the event in CLAS and hit in the tagger counter, 
DOCA cuts on the A and the S~ candidates, detached 
vertex cuts for the reconstructed A, and particle identi- 
fication cuts on /3 versus p. The choice of these analysis 
parameters did not affect the final result for the upper 
limit of the cross sections, and we therefore considered 
their uncertainties to be negligible. 

The uncertainty in the overall normalization of the ex- 
periment is driven by the relatively large trigger efficiency 
corrections that are required, as described previously. We 
assign an uncertainty of ±15% to the absolute normaliza- 
tion, which is estimated by comparing our determination 
of A++ production to other measurements, as discussed 
in the previous section. 

The dominant contribution to the photoproduction 
cross section systematic uncertainty comes from the 
model dependence of the estimated acceptance, and we 
take the RMS of the model calculations to obtain a ±21% 
relative uncertainty. 

Table IIVI summarizes the dominant contributions to 
the systematic uncertainties. We used the sum in quadra- 
ture of 26% to determine the upper limits for the exper- 
imental cross sections described in the next section. 



UPPER LIMITS 

In order to determine the upper limits for the cross 
section for a possible peak, we scanned each 6 MeV-wide 
bin in the sideband-subtracted mass spectrum in Fig. [HI 
considering the center of each bin as the mean value of a 
Gaussian distribution with a fixed width oi ac — 1 MeV, 
which represents our experimental resolution for a poten- 
tial $ (1860), as estimated by the detector simulation. 
Then the points in the neighborhood of the bin are fitted 
to a Gaussian peak plus a polynomial. The purpose of 
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SUMMARY 
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FIG. 9: [Color online] Upper limits for the photoproduction 
cross section with subsequent decay ^ — > H~7r~ in 20 MeV 
mass windows at 90% CL versus the invariant mass of 3~tt~ . 
The red dashed vertical line marks the position of the en- 
hancement reported by NA49 . 



this x^-fit is to provide us with an estimate of the back- 
ground under the possible peak and its uncertainty. The 
total number of signal events is calculated as the excess 
of the observed events over the fitted background, both 
integrated within a window of 20 MeV around the cen- 
ter of each bin. The red dashed curve in Fig. [S] shows a 
particular example of the fit for a mass bin centered at 
Afc-^- = 1.813 GeV. 

In order to obtain an upper limit on the photoproduc- 
tion cross section, we developed a procedure based on the 
method described in Ref. [40] ■ This procedure allows us 
to factor in the uncertainties in the background extrac- 
tion and the acceptance correction into the determination 
of the upper limits at a given confidence level (CL) using 
the sideband-subtracted method. We also performed a 
cross check of our method with an approach for estimat- 
ing the upper limits by Smith 4l| based on the construc- 
tion prescription of the Feldman- Cousins method [4^ . 
which properly takes into account the systematic uncer- 
tainties when constructing the confidence belts. For a 
comparison, we assumed a mass-independent acceptance 
of 0.4% and a mass-independent relative acceptance un- 
certainty of 30%; the agreement between the two meth- 
ods is very good. A more detailed description of our 
method to determine the upper limits can be found in 
Appendix A. 

The upper limit of the photoproduction cross section 
at 90% confidence level for the process — >■ $ X — >■ 
S~7r~Ar versus invariant mass of the S~7r~ system is 
shown in Fig. [HI In the mass range near M = 1.862 GeV, 
where the NA49 collaboration observed an enhancement, 
we obtain a 90% CL upper limit of ^ 0.7 nb. 



In conclusion, we have conducted for the first time a 
search for the $(1860) pentaquark state in real pho- 
toproduction within the incident photon energy range 
4.5 GeV < Ej < 5.5 GeV. We do not observe any 
statistically significant enhancement near invariant mass 
M = 1.862 GeV. The upper limit at 90% confidence 
for the photoproduction cross section of the reaction 
'yd — > ^^^X multiplied by the branching ratio for 
$ — > S^TT^ is determined as a function of the in- 
variant mass of S~7r~, using a method similar to the 
one described in Ref. [40|. The upper limit for the cross 
sections for $ photoproduction with subsequent de- 
cay $ — > S^TT^ for 20 MeV mass windows is less 
than 3 nb in the S^tt^ mass range between 1.6 GeV 
and 1.9 GeV. The upper limit is less than 1.5 nb for the 
masses from 1.9 GeV to 2.2 GeV. The upper limit for 
the cross section averaged within a narrow mass range 
of 1.80 GeV < Mh^ < 1.92 GeV is ~ 0.7 nb. This 
is approximately a factor of three improvement over the 
previously estimated upper limit in small-angle electro- 
production by the HERMES collaboration of 2 nb 0]. 
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Appendix: Determination of upper limits 

The number of events in the 20 MeV mass windows 
from the sideband-subtracted spectrum in Fig. [5] are 
distributed according to a Gaussian distribution with a 
width determined by the statistical uncertainty obtained 
during the sideband subtraction. Therefore, in each mass 
window we model the excess of the events above the back- 
ground, or signal events, according to a Gaussian distri- 
bution, with a mean fJ. > limited by the condition that 
the cross section cannot be negative even if the number 
of observed events is less than the expected background. 

Systematic uncertainties are included into the calcu- 
lation by assuming that the measured acceptance and 
the number of background events are random variables 
distributed according to the normal distribution: 



P{x,bm,em I fJ;b,e) 



e ^ 



2tT(7x 

1 



(3) 



2'Kae 



where /i is the expectation of the number of signal events, 
e is the acceptance factor for the signal, b is the expec- 
tation value of the number of background events in the 
20 MeV mass window, and x is the observed number of 
events in the same window for each fit. The experimental 
statistical uncertainty is denoted by ax- The estimated 
number of background events bm is determined from the 
polynomial fits to the background, as shown in Fig.[6l and 
its uncertainty, df,, is taken from errors returned from the 
fit. The estimated signal acceptance Cm is determined by 
Monte Carlo and shown in Fig. [H For the uncertainty in 
the value of e,„, Ue, we use the systematic uncertainty of 
26% (Table [r7|. P(a;,5™,e™ | M,5,e) is the probability, 
under assumption of our model, to observe the values x, 
bm and Similar to the Rolke method in Ref. [4o| we 
use a profile likelihood to estimate the confidence level 
(CL). The logarithm of the profile likelihood is defined 
as the logarithm of the ratio 



C = 



-21nA(^t) 

sup{P(a;,5„i,e,; 



-2 In 



fit,b, e);b,e} 



^sup{P(2;, 



fj.,b, e);fj.,b, e} 



(4) 



where A is the profile likelihood, and /it is the hypothesis 
value being tested. The supremum, or the least upper 
bound, sup{P{x, bm,Gm I /it, b, e); ^, b, e} in the denomi- 
nator under the logarithm in Eq.[4]is taken over all values 
of {ii,b,e), and is located at (6 = 6m, e = Cm,// = ^^). 
The least upper bound sup{P(x, bm, e™ | /i, b, e); b, e} in 
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FIG. 10: [Color online] Illustration of determining the upper 
and lower limits for CL = 90%. The a;-axis is the expectation 
value for the measured quantity /x. The red solid curve is the 
logarithm of the profile likelihood C versus , and the green 
dashed line is ~ 2.705 corresponding to CL — 90%. 



the numerator is taken only over the background and ef- 
ficiency h and e. In order to determine the location of 
this supremum we find the zero-crossings of both par- 
tial derivatives of the likelihood P in Eq. HI ^ = 



and ^ = 0. In the limiting case when the background is 
known to be zero and the efficiency is 100%, Eq.|3]simply 
becomes 



C 



(5) 



The log-likelihood distribution in Eq. 2] is approxi- 
mated by the x^-distribution (x^ « £) with the appro- 
priate number of degrees of freedom. In this case, there is 
only one degree of freedom. To find the values of /xt corre- 
sponding to a certain confidence level, one first finds the 

value corresponding to that CL for the x^-distribution 
with a single degree of freedom. The solutions for fit 



are the values where C differs from its minimum by that 
amount of x^- The solution that is less than the most 
likely value of /i is the lower limit, while the larger so- 
lution is the upper limit. After substituting x^ for ^ 
in Eq. (U) and solving the quadratic equation, the two 
solutions for fit for a given x^ can be found as follows: 



2 2 
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2 2 2 



(6) 



4 2 2 



2 27,2 



4 2 2\ 

X cre<^b) 



(22 
[X CT, 



Fig. [TU] illustrates how the upper and lower limits at 
90% CL are found. The red curve is the log-likelihood C 
with a minimum at around fXt ~ 3300, which is the most 
likely value for the illustrated example. The probability 
of having 90% of the trials within a certain range is re- 
alized for x^ = 2.705, which is represented by the green 
line in Fig. (TU] The intersection points of these curves 
give us the upper and lower limits at 90% confidence 
level. In certain cases, for instance when the estimated 
efficiency is very low, and the uncertainty for it is rela- 
tively large, no upper or lower limits can be found. This 
happens when the logarithm of the profile likelihood C 
does not behave like a parabola, and therefore is not well 
approximated by the x^-distribution. 

In cases where the most likely /if is negative, we take 
the lower limit to be 0. If the most likely value fit and the 
upper limit are negative, we increment the number of ob- 
served events x by one unit until we get the first positive 
value. Such an ad hoc adjustment can cause the coverage 
probability, defined as the probability that a true value 
of the cross section for a process is less than the corre- 
sponding upper limit obtained by this method, to differ 
from the desired confidence level of 90%. Therefore, one 
needs to check that the results obtained by this procedure 
indeed provide the desired confidence level. Our Monte- 
Carlo tests showed that the coverage probability of this 
method is within 5% of the nominal confidence level. 



